Stochastic Resonance: from climate to biology 
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In this paper I will review some basic aspects of the mechanism of stochastic resonance. Stochastic 
resonance was first introduced as a possible mechanism to explain long term climatic variation. Since 
then, there have been many applications of stochastic resonance in physical and biological systems. 
I will show that in complex system, stochastic resonance can substantially change as a function of 
the "system complexity" . Also, I will briefly mention how to apply stochastic resonance for the case 
of Brownian motors. 

PACS numbers: 



1. INTRODUCTION 

The first numerical simulation providing strong evi- 
dence of Stochastic Resonance was performed by An- 
gelo Vulpiani and myself in a rather exciting night on 
February 1980. Together with Alfonso Sutera and Gior- 
gio Parisi, we were trying, at that time, to understand 
whether a relatively small periodic forcing can be ampli- 
fied by internal non linear stochastic dynamics, leading 
to a possible understanding of the 100 Kyear cycle ob- 
served in climate records. The first version of our paper 
was not accepted for publication in Journal of Geophy- 
isical Research and in Journal of Atmospheric Science. 
Eventually the paper was published in Tellus in 1982 
together with a similar paper by Nicolis and Nicolis 
who, independently, were proposing the same mechanism 
for climatic change. Short after that february night in 
1980, we were able to provide a quite general understand- 
ing how the mechanism works and how to generalize it 
for chaotic systems. Both theory and generalization ap- 
peared as a letter in J. Phys A in 1981 i.e. one year 
before the Tellus paper! 

We were, and still we are, convinced that Stochastic 
Resonance is a rather new and conceptually important 
phenomenon in science. For almost a decade, only few 
colleagues were shearing with us the same feeling. Even- 
tually, after the paper of Moss et.al. [4] and some new 
theoretical work [5[ stochastic resonance became a fash- 
ionable and interesting research topic in many different 
scientific areas, see [6| for a review. 

Resonance is a kind of magic word in physics and there 
is no surprise that a new resonance mechanism can excite 
the scientific community. In our case, the name stochas- 
tic resonance was introduce because of a short discussion 
with J.Imbrie on march 1980. Just few days after our first 
numerical simulation, I participated to a climate meeting 
in Erice where I gave a short talk on our results. J.Imbrie, 
one of the most famous scientist working on paleoclimate, 
asked whether what we found was somehow similar to a 
resonance and my answer was: "not exactly! It is a kind 
of stochastic resonance!" . This is how the name was in- 



troduced in literature. In some sense, we can think of a 
"resonance" as follows. In the standard resonance mecha- 
nism, think for instance of a damped harmonic oscillator, 
the amplification of the external forcing can be related, 
mathematically, to a singularity in the complex plane 
of the green function of the problem. The real part of 
the singularity is, of course, the resonance frequency. In 
chaotic or stochastic systems, there are singularities in 
the complex plane but on the imaginary axis. The mech- 
anism of stochastic resonance provides a way to shift the 
singularity on the complex plane with a non zero real 
part. This rather cumbersome view of stochastic reso- 
nance may be not entirely clear but it tries to justify 
why the word "resonance" can still be used beside any 
historical reason. 

In this paper, I will review the main idea on stochastic 
resonance with an emphasis on climate research, where it 
was originally proposed. I will also show how stochastic 
resonance is still an interesting subject to work on and 
in particular I will discuss some new findings in simple 
and complex systems. Finally, I will mention how mech- 
anisms related to stochastic resonance are of interest in 
biology and in the wide research field of brownian motors. 



2. THE MECHANISM OF STOCHASTIC 
RESONANCE IN CLIMATE THEORY 

Climate is one of the most complex system in nature 
and it is a major scientific challenge to understand the ba- 
sic mechanisms leading to climate changes. In the past, 
climate shows quite remarkable changes over all times 
scales. In particular, one of the most striking aspects 
of past climate changes is the so called lOOky cycles ob- 
served in paleoclimatic records. The cycle itself is not 
at all periodic showing sudden warming phases followed 
by gentle decreases of temperature. This saw-tooth be- 
haviour is in phase with the so called Milankovitch cycle, 
i.e. to the change of global radiation of the Sun due to as- 
tronomical change in Earth orbit. This is the only global 
(i.e. averaged) change experienced in the incoming radi- 
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FIG. 1: The energy balance climate model. In the upper layer 
we show the basic variables of the model: the incoming radi- 
ation Rin and the out coming radiation Rout and the infrared 
emission. In the lower panel we show the albedo as a function 
of temperature Rin — Rout = (1 — a{T))Rin- The thick line 
represents the infrared emission. 

ation. Is that possible to relate the astronomical forcing 
to the observed climate change? This is the basic ques- 
tion which puzzles scientific research since many decades. 
A simple computation shows that it is not easy to guess 
the correct answer. 

Let us consider the simplest possible climate model. 
We consider the averaged Earth temperature T as the 
basic variable we want to describe, the precise meaning 
of average is of marginal interest for the time being. We 
know how much radiation is incoming, let us say Rin- 
The out coming radiation Rout depends on the earth's 
temperature T (infrared emission) and reflection from 
the surface, which can be written as aRin- Assuming 
Ce the earth's thermal inertia, we can write the energy 
balance model as: 

dT 

Ce-^ — Rin ~ Rout (1) 

where Rout = aRin-\- Ej^ Ej being the infrared emission. 

What we might call Earth'climate is just the station- 
ary solution of equation ([T]). While Rin is independent 
on T, both a and Ej should be dependent on the Earth's 
climate which, in our simple case, is described by the av- 
eraged temperature T. In a first approximation, we can 
think that Ej is linearly dependent on T (a more com- 
plex behaviour does not substantially change the results). 
The so called albedo a is also a T dependent quantity. 
More specifically, we expect a to be quite large for ice 
regions. An appropriate way to think of the function 
a(T) is shown in figure ([T]). Following the cartoon shown 
in the same figure, it is clear that there exist more than 
one stationary solution of equation ([T|), corresponding to 
three different possible Earth climates. The coldest of 



the possible climates is referred to as "ice cover earth" 
climate, i.e. the case where almost all incoming radia- 
tion is reflected back in the space. A careful estimate 
of the albedo effect and the infrared radiation, obtained 
by current observations, shows that the warmest possible 
climate is close to the current Earth climate (we discard 
the intermediate climate simply because it is dynamically 
unstable). Thus equation ([T]) gives us a quite amazing re- 
sult as it explains our present climate state. Note that 
an increase of CO2 in the atmosphere is equivalent, in 
our simple model, to a change in Ej represented by the 
thick line in the lower panel of figure ([T]). 

Our simple model tells us that there is an interesting 
feature in climate dynamics represented by the instabil- 
ity of the intermediate climate state. What is the phys- 
ical origin of this instability? The answer is related to 
the complex mechanism of cloud formation. Clouds by 
themselves reflect incoming radiation. When the tem- 
perature increases, the albedo start to increases as well 
because more clouds are expected to form. Eventually, 
for large enough temperature, the amount of cloud for- 
mation is balanced by precipitation leading to a constant 
value of a. For the time being, let us concentrate on the 
present (stable) climate. The orbital forcing, discovered 
by Milankovitch, is a small amplitude modulation of Rin 
in time, i.e. 

Rin{t)=Rf^^Acosiujt) (2) 

Let us consider the effect of the Milankovitch forcing in 
equation ([T]). If T is always quite close to the stationary 
(present) climate Tq, we can estimate ST = T — Tq as: 

^ = --ST ^ Acos{cot) (3) 
at T 

where 1/r is the relaxation time of the present climate. 
More precisely r can be estimated by the equation: 

r Ce dT 

After an appropriate estimate of A, one can easily show 
using (|3]) that the effect of the orbital forcing on the 
present climate is of order At ~ 0.5i^, much smaller 
than the lOi^ observed in paleoclimatic records. Thus, 
it seems that orbital forcing can hardly explain the 
glacial/interglacial transition experienced by the Earth 
climate. 

In order to make progress we need to supply our simple 
model by new features. An important feature is to as- 
sume that there exist more than one stable climate state 
close to our present value Tq. More precisely, we shall 
assume that the albedo feedback mechanim (described 
in terms of clouds and ice) can eventually lead to a more 
complex behavior near our present climate leading to two 
stable climate states Ti and Tq separated by about 1{)K 
see figure ([2]). This is a rather ad hoc assumption which 
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is, at any rate, consistent with observations on albedo. 
Still, this feature does not really change our conclusion 
since the orbital forcing is too small to provide a suitable 
mechanism to produce transitions between Tq and Ti. 

The situation can improve dramatically if we consider 
the following three non trivial statements: 

• Many complex systems can be described by means 
of slow variables and fast variables, even if there is 
no explicit time scale separation between the two 
kind of variables. Fast variables may be considered 
as "noise" acting on the dynamics of slow variables. 

• For Earth climate, the averaged temperature in ([T|) 
should be considered as a slow variable with respect 
to the fast variables due to weather variability. 

• When there exist multiple equilibriums, the noise 
can introduce a long (random) time scale to switch 
from one equilibrium to another. 

The three statements are no trivial and should be ex- 
plained carefully. For many decades, the splitting among 
slow and fast variables have been considered a well known 
feature of molecular motion. At the atomic scale, parti- 
cles continuously experience short time interactions with 
other atoms while the hydrodynamic behavior (the slow 
variables) is described by a suitable space average over 
the molecular chaos. Looking at the fluctuating proper- 
ties of a turbulent flow, there exists a well deflned scale 
separation (both in time and in space) from molecular 
motion and hydrodynamic flow. Nevertheless, even in a 
turbulent flow, i.e. disregarding any molecular motion, 
we can still distinguish large and small scale fluctuations. 
In this case, however, there is no spectral gap identify- 
ing any scale separation among scales. Therefore, one 
may be tempted to state that, in this case, the concepts 
of fast and slow variables are poorly defined. It turns 
out that such a statement is too limited. There exist 
experimental systems where the large scale motion can 
be thought as the slow variable superimposed to a noise 
(turbulent) small scale background acting on the system. 
Recently [7], it has been shown that the above picture 
can explain qualitatively and quantitatively the numeri- 
cal results observed in a simplified model of turbulence. 
Thus, in most complex systems the separation of large 
scale, slow variables and small scales, fast variables does 
correctly picture the dynamics in a self consistent way. 

If we apply the above discussion on the Earth climate 
(second item), then on the time scale corresponding to 
the orbital forcing, the day by day weather fluctuations 
should be considered as noise (fast) variables acting on 
the climate system. This is not just a way of thinking. 
Disregarding the effect of noise is a major limitation for a 
correct description of the physical properties of climate. 
It took a long time, and it is still on going, for the sci- 
entiflc community to accept that noise is not a measure 



of our "ignorance" but just a physical feature of most 
complex systems. 

The third item is more delicate and needs to be ex- 
plained in details. Let us assume that our climate system 
(although the discussion is true for any system) shows 
two stable equilibriums whose temperature difference is 
2AT. In this case, it is well known that the effect of the 
noise can induce transitions between the equilibria. The 
characteristic time tl for the transition can be estimated 
as: 



Texp{AT^/aT) 



(4) 



where as in equation ([3|) r is deflned as the (fast) relax- 
ation time over one of the stable equilibriium and cr is the 
variance of the noise (weather fluctuations) acting on the 
system. Thus our year by year temperature is assumed 
to be decomposed in the climatic component T and a fast 
noise component. In equation (|4]), tl should be consid- 
ered an average time to switch from one equilibrium to 
another. In other words, the transition time between the 
two equilibriums T^q is a random variable with average 
value Actually, under some quite generic conditions, 
we can predict the probability distribution of Teq which 
is 



exp{-Teq/TL) 



TL 



(5) 



Let me summarize the new feature we have introduced in 
the simple energy balance model ([T]). We have assumed a 
new form of the albedo which is providing three different 
climate states around To, let us call them Ti (stable), 
Ti (unstable) and Tq. The shape of the new albedo is 
shown in flgure ([2]). We assume that Tq — Ti = 1{)K and 
Tq — Ti = Ti — Tl. Next we have introduced an external 
noise which takes into account the dynamics of the fast 
variables, whatever they are, not represented by the slow 
component T. We can rescale T near Ti and simplify the 
equations by using T = Ti ^ ATX with AT = 5^i^: 



dX_ 
~dt 



(6) 



where r]{t) is a gaussian random variable with unit vari- 
ance {{r]{t)r]{t')) = S{t — t')) and a is the noise strength. 
Notice that we have also rescaled the time unit. In order 
to do some computation what we still need is the strength 
of the noise. We can obtain some information by the 
following trick. In our present climate (i.e. To) we have 
observational data which tells us about average tempera- 
ture fluctuation. Also, we know from General Circulation 
Models, which is the time correlation for temperature 
fluctuations. Thus, we can estimate the strength of the 
noise by linearizing our climate model around Tq and us- 
ing the fluctuation-dissipation theorem, i.e. temperature 
variance = noise strength / (2 relaxation time). Once we 
know the noise strength cr we can compute the average 
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FIG. 2: The albedo feedback for an energy balance model 
with three climate states near To. The model is conceptually 
similar to the one described in figure ([!]) with a change in the 
albedo. 

switching time from Tq to T\. It turns out that, by using 
(j4|) the average transition time between the two stable 
states of our climate model is close to 50, 000 years! This 
is a rather surprising and interesting result because it 
tells us that there exists a rather long time scale (50, 000 
years) introduced by the nonlinearity in the model and 
the noise, to be compared against the quite short time 
scale of the model, namely the deterministic relaxation 
time towards one of the two stable equilibriums which 
is of the order of 10 years. However, this is not enough 
because the characteristic time of 50, 000 years is a ran- 
dom time with an exponential probability distribution. 
Therefore nonlinearity and noise are not enough to get a 
periodic behavior of the temperature. 

In order to get a periodic behavior, the idea is to in- 
troduce the Milankovitch effect in the system. After a 
simple computation one gets the following results: 

^ = X - + Asin{yjt) + ^/Gr](i) (7) 

An appropriate rescaling in terms of physical quantities 
gives us the value oi A ^ 0.11 for the Milankovitch terms 
andcj = 27r/10^xl/10. Without the noise, the value of X 
changes periodically in time with an amplitude order A. 
Going back to the temperature, this implies that, with- 
out noise, we have a periodic behavior of temperature 
with amplitude AAT ~ 0.5K. The situation changes 
completely when we consider the effect of the noise. A 
correct understanding of the noise effect can be achieved 
by considering the simplified equation: 

'^ = x-x^^A^n{t) (8) 

i.e. without any time dependency on the external forcing. 
It is crucial to remind that for any stochastic differential 
equation of the form: 



FIG. 3: Numerical simulation of equation (|7|) in three different 
cases. In upper panel, we show the solution with the noise and 
A — ^. In the middle panel we show the perturbing forcing 
Acos(2i^ujt) while in the lower panel we show the solution 
with A 7^ 0. The effect of the small periodic forcing is to 
synchronized the random switching from one climate states 
to the other, i.e. to get a stochastic resonance. 

the "equilibrium" probability distribution P(X) of X is 
given by: 

P(X) = Z^exv{-2V{X) /(j) (10) 

where Zn is a normalization factor. In our case V^(X) = 
-1/2X2 ^ 1/4x4 _ ([HI) and we obtain from ^ 

P(X; A) = P(X; A = {))exp(2AX/ a) (11) 

The probability distribution when A ^ changes dra- 
matically the probability distribution when A = 0. Even 
for small A the exponential factor exp{2A/a) can be 
rather large, i.e. order 30 in our case. The consequence of 
this simple calculation is that the probability to be near 
X ~ 1 is 30 times larger than to be in X = —1. Now it is 
quite clear that if A is slowly varying as cosiuot) in time 
than the probability distribution is peaked periodically 
in ±1 with almost the same period of 2i\ juj. This is the 
very essence of stochastic resonance. 

In figure (|3]) we show the numerical simulation of ([7j) 
and in particular the upper panel is X(t) for A = 0, the 
pnnel in the center shows the behaviour of A(i) while 
the lower panel is the solution of ([7j). As we can see, the 
effect of the small periodic perturbation is to synchronize 
the "random" behavior of X and to produce an almost 
periodic output. In some sense, stochastic resonance is a 
counter intuitive phenomenon because without noise the 
system shows a small amplitude modulation around one 
of the two stable steady states ±1, while adding the noise 
we obtain a large amplitude effect with the same period. 



(9) 



One has to be a little bit carefully in understanding the 
possible results of ([7j). As we said, the modulation of A in 
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time should be slow enough in order to get a stochastic 
resonance. We can be more quantitative using the fol- 
lowing simple argument. Let us go back to equation (|8|). 
We can compute the probability distribution of the ran- 
dom time to go from +1 to —1 and we already know that 
the probability distribution is exponential. Let ri(+l) 
the average exit time, then the variance of the random 
switching time is also ri(+l). Finally we know how to 
relate ti(+1) to A, namely by using the expression: 



exp{-{AVo 
a 



A)) 



(12) 



where in our case AVq = 1/4. Thus the requirement we 
need is that ri(+l) <C 27r/a; and ri(— 1) ^ 27t/u; where 
Ti(— 1) is the average exit time starting with the initial 
condition X = —1. The two requirements implies 



27r 
27r 



exp(-(Al/o - ^)) < 



(13) 
(14) 



The meaning of the above inequalities is represented in 
the cartoon of figure (|4]). The time is in the counter- 
clockwise starting with the left top panel. We start at 
X = 1 and as time goes by, the system at t = T/2 has 
a very high probability to jump in X = —1 if ([Tlj) is 
satisfied. Then at t = T we have the same effect for the 
transition X = —1 ^ X = 1. The two conditions (|13|14p 
tell us that if the noise is too small then we will never 
see anything and, on the other hand, if the noise is too 
large transitions will happen regardless of the external 
periodic forcing A. Thus, there is a range of a (i.e. the 
resonant range) where we obtain periodic transitions in 
phase with the external forcing. 

One can highlight the resonance effect also by using 
another kind of variable. Let X{t) the solution of ([7j) 
and let us denote Fx{i^) its Fourier transform. Then 
represents the power spectrum of X. We should 
expect that \Fx{i^ = cc;)^ is a function of the noise with 
some maximum value for the "resonant noise" ctr. This 
is indeed the case as shown in figure ([5]), where \Fx{i^ = 
cc;)p is plotted against the a, while in the insert we show 
\Fx{iy)\^ for cr = (JR. 

Let me know go back to the climate. As we have seen, 
the effect of a small periodic forcing in our "simplified" 
climate model is to produce a periodic glacial to inter- 
glacial transition in phase with the Milankovitch forcing. 
This is a quite good new. However, the signal shown in 
figure ([3]) poorly compares against the observed behav- 
ior of the "proxy" earth temperature. The latter shows 
a kind of saw tooth behavior with strong asymmetries 
between cold and warm period. Thus it seems that our 
theory cannot represent the real climate. Such a nega- 
tive statement, however, is somehow misleading simply 
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FIG. 4: A cartoon to explain the mechanism of stochastic 
resonance. The period of the oscillation is T. At time t = 
the climate states has low probability to jump from +1 to —1. 
At time t = T/2 the situation is reversed: the probability no 
to jump is extremely small. At time t = T we start a new 
cycle. 



because the model so far used, namely the energy bal- 
ance model (j6j), includes only the simplest possible feed- 
back due to radiation. More sophisticated effects, like 
the temperature-precipitation effect, can change the sit- 
uation. Based on the above discussion we can make the 
following statement: 

• stochastic resonance is able to show that the exter- 
nal "weak" forcing due to Milankovitch cycle can 
be amplified by internal non linear dynamics of the 
climate; 

• climate dynamics shows a long time scale behavior 
due to non linear feedback; 

• noise can develop long time scale behavior if there 
exists multiple equilibriums; 

The most important feature coming out from our analysis 
is that it is crucial to understand the non linear interac- 
tion between fast a and slow time variables in climate 
models. Time scales of order 1000 years or more, do no 
rule out the importance of the dynamics of fast variables 
which usually characterize atmospheric and/or oceanic 
circulations'. This, I think, is an important feature even 
in more sophisticated and/or realistic climate models. 

The mechanism of stochastic resonance has been re- 
cently invoked to explain the observed climatic variation 
on a shorter time scale with respect to the Milankovitch 
forcing. During the last 100, 000 years, the earth climate 
exhibits abrupt changes from relatively cold to warm 
states, different to those observed in the 100, 000 year 
cycles. According to Alley et.al. [1], the climate changes 
occurred on time intervals which are multiple of 1,500 
years, i.e. the time to switch from one climate state to the 
other is either 1,500 or 3,000 or 4,500 and so on. This 
rather peculiar behavior calls for an explanation which. 
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not surprisingly, can be given in terms of stochastic res- 
onance. As we previously discussed for equation ([7j), the 
effect of the periodic forcing is to decrease and increase 
(depending on the time t) the probability of a switching 
from a climate state to the other. With reference to fig- 
ure (|4]), if at time t = T/2 the system does not switch, 
then it takes a full period T for the possibility to make 
a jump. This is certainly the case if the noise intensity 
is smaller than ctr. It follows that the probability of the 
exit time is not peaked around T/2 but it is quantized 
with maxima in (2n + l)T/2, n = 1,2, .... To show that 
our analysis is correct, in figure (|6]) we plot the proba- 
bility distribution of the time to switch for a = O.SaR 
(symbols) and for a = cjr (continuous lines). The quan- 
tization effect is quite clear as predicted by our simple 
analysis. 

It appears that the observation of Alley et.al. may be 
explained by the mechanism of stochastic resonance as it 
was pointed out by the same authors. Actually, a theo- 
retical analysis of a simplified ocean/atmospheric model, 
due to Ganopolski and Rahmstorft Q shows that the 
effect of a periodic forcing in freshwater input over the 
North Atlantic Ocean combined with a suitable stochas- 
tic forcing, produces a stochastic resonance in cold/warm 
climate changes as observed in real data. In the Ganopol- 
ski and Rahmstorft model, the two climate states, be- 
tween which one has climate transitions, represent two 
different thermoaline circulation in the Norht Atlantic 
Ocean, which are responsable for sea ice growth and de- 
stroy. 

It is unclear, at this stage, whether stochastic reso- 
nance is or is not the possible explanation of the Alley 
et.al. data and more research is needed to be done for 
a reasonable assessment. At any rate, it is important to 
remark that abrupt climate changes as those observed 
during the Milankovitch cycle or over the last 100, 000 
years are hard to be explained without the idea of transi- 
tion between multiple equilibriums triggered by the noise, 
which is the basic feature of stochastic resonance. 



3. STOCHASTIC RESONANCE IN COMPLEX 
SYSTEMS 

In the last ten to fifteen years the mechanism of 
stochastic resonance has been widely applied to a num- 
ber of different physical and biological systems. Among 
them, the most striking applications concern neural sys- 
tems pioneered by Moss and collaborator In most 
cases, experimental results show a remarkable agreement 
with respect to the qualitative and quantitative picture 
represented in figure ([6]). Also, in many applications one 
does not know the "equation of motions" for the system 
and, actually, periodic perturbation leading to a stochas- 
tic resonance are used as a "measure" of non trivial co- 




FIG. 5: The figure shows the Fourier amplitude |Fx(z^ = 
for the solution of equation d?]) for different values of the noise 
amplitude a. In the insert we show \Fx{i^)\^ for a = cfr^ 
i.e. the optimal noise for which the Fourier amplitude at zy = cj 
is maximum. 
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FIG. 6: Probability distribution of the random switching time 
r for two different noise amplitude obtained as the solution 
of equation (|7|). The amplitude and the period of the forcing 
is the same. Line refers to gr (see figure (|5|) while symbols 
refer to a = O.Scfr. In the insert we show the same figure 
in log-lin plot. It is evident that for a smaller than gr the 
probability distribution of r is "quantized" at integer values 
of (2n + l)T/2. 

operative phenomenon in the system under study (this is 
certainly the case for most neurophysiological systems). 

When dealing with complex systems, one can find non 
trivial behavior of the stochastic resonance. In this sec- 
tion I show that this is the case for a network system 
composed by some "node' ^/^i, i = 1, ...A^, see figure ([7]). 

Just to simplify the discussion, I assume that the dy- 
namics of each node is controlled by the same equation 
and that the connectivity of the system is described by 
the matrix Lij'. 

^ m^l)i - gi/jf + Liji/jj + y^rji (15) 

where m > and g are real constant. Also, I will assume 
that ^iLiji/jj = which clearly can be done without 





FIG. 7: A cartoon of our "complex" system. Each node sat- 
isfies equation ()15p . Lines describes the connectivity matrix 



lack of generality, and finally, that Lijipi^pj > 0, i.e. 
the connectivity matrix does not introduce any "insta- 
bility" in the dynamics of the system. The noise rji is 
(5-correlated in time and {r]ir]j) = Sij. 

The complexity of our problem, so to speak, is intro- 
duced by the matrix Lij. We want now to study the be- 
havior of (p!5]) when an external periodic forcing is added 
to the system, i.e. 



dt 



= rmpi - gijjf + Lijijjj + ^rji + AiCos{2ujt) (16) 



In order to simplify our work, let us focus on the "aver- 
age" variable ^ = ^-ipi. One should expect that 
for long enough period T = 27r/cj and a suitable noise 
(jR{m^g) a stochastic resonance can be observed. Note 
that I have defined ctr as an explicit function of the vari- 
ables m and g. By averaging (fT5]) we obtain 



~dt 



(17) 



where ^ = N'^ and (...) = N'^ ^.i I as- 

sume that e is independent of A/', i.e. a is chosen in such 
a way that e is a fixed quantity. The difficulty is to com- 
pute the term (V^f). For this purpose, we define such 
that tpi = <l> + i.e. (f)i are "deviation" of ipi from <l>. 
Then we have (tl^f) = <l>^ + 3g^{(j)f). This expression is 
correct as far as we can neglect the term ((/)f ) which, in 
most cases, is a good first approximation. Putting all 
together, we obtain: 



~dt 



(m - ?>g{(fl))^ - g^^ ^ yfer] 



(18) 



As we can see, the effect of complexity, i.e. the con- 
nectivity matrix L^j, introduces a change in the linear 
term which now become a time dependent function . It 



FIG. 8: Numerical simulation of equation ()15|) with the same 
parameters m,g,a and two different topology (shown in the 
figure): one dimensional lattice (upper panel), two dimen- 
sional lattice (lower panel). The amplitude and the period 
of the periodic forcing is the same for both cases. While the 
two dimensional lattice shows stochastic resonance, the same 
is not true for the one dimensional case. 



can happen that (m — 3g{(j)'^)) < and if this is the 
case, transition between the two states (whatever they 
are) occur with a mechanism completely different with 
respect to what we discussed in the previous section. Is 
this the case? Everything depends on the quantity (0f ) 
and therefore on the connectivity matrix Lij. 

In order to make progress, let me first discuss what 
are the "statistical stable" equilibriums of (p!8|) . Let us 
define 



(19) 



We should expect that the statistical stable equilibri- 
ums are ±^0- Then fluctuations around the equilibrium, 
which in first approximation we can describe as 0^ , satisfy 
the equation: 



dt 



L, 



(20) 



where 



a={m- ig{cl>1)) - ig^l = -2(m - 3g{<l>^)) (21) 

Using ([2Q|) we can now estimate Let us define — 

the eigenvalues of Lij. Then, by using a well known 
result in the theory of stochastic differential equation, 
we obtain: 



a + An 



(22) 



In the limit of large A/", we can define the density of states 
p{X) with the approximation: 



d\p(\) 



(23) 



Putting everything together, we finahy have: 



1 + Ar 



2 ./ 2(m 



P(A) 



(24) 

Equation ([24|) is a non hnear equation relating the value 
of to the density of states p(A),i.e. to the connec- 
tivity matrix Lij. Once we have we can comput <l>o. 
The tricky part of our problem is that we have computed 
((/)|) when (^4) = 0, i.e. with no external forcing. If we 
now have an external forcing, all our computation for 
and a should change taking into account A. One can 
compute perturbatively the effect in power of A and take 
the first order for small A. The computation are done for 
a special case in [10]. The final results is that everything 
goes as in the theory discussed in section 2 but with a 
renormalized value of A, i.e. 



^(1 



l-2Dg 



where 



D 



d\- 



P(A) 



(25) 



(26) 



(ao 

and Qfo corresponds ([2T]) for A = 0. For any practical 
purpose our result means the following. Let us imagine 
a network described by equation (p!5]) . Then the effect 
of a periodic forcing with amplitude A on the average 
$ is equivalent to a one dimensional problem (similar to 
(|6|) with a renormalized amplitude ([25]) . The connectiv- 
ity matrix fixes the value of the renormalization by (|26|) . 
Thus, depending on L^j, the effect of stochastic resonance 
can be enhanced or depressed. In figure ([8]) we demon- 
strate our result in a simple case. We consider the same 
system (i.e. same value of m^g^A and <j) for connectiv- 
ity matrix topologically equivalent to a one dimensional 
lattice (upper panel), and a bidimensional lattice. In the 
latter case the system shows stochastic resonance while 
the same is not true for the former case. 

One main conclusion that we can outline from our dis- 
cussion is particular relevant for climate theory. The 
physical meaning effect of an external forcing can be dras- 
tically different depending on the feedback in the system 
(in our case the matrix Lij). There are cases, as we have 
seen, for which a small forcing or a small noise can trig- 
ger transitions and a naive computation of the relevant 
quantities (i.e. neglecting the effect discussed in this sec- 
tion) can lead us to wrong results. I want to argue that 
the above conclusions is relevant for other physical and 
biological systems. 



4. BROWNIAN MOTORS AND STOCHASTIC 
RESONANCE 

Up to now we have discussed the case of stochastic res- 
onance in a more or less traditional fashion, i.e. when 




cold 



FIG. 9: A simple model of Brownian motor. Going from the 
top to the bottom of when the temperature is low (cold) a 
particle moving in the ratchet potential is trapped in one of 
the minima. Rising the temperature, the particle becomes 
almost free to move and, because its initial condition (mini- 
mum) is close the maximum on its right, the particle has more 
probability to go towards the right minimum. If the temper- 
ature is now decreasing the particle has more probability to 
be trapped in the right minimum with respect to be trapped 
in the minimum on the left 



multiple equilibriums exist. We want now to under- 
stand whether there are other physical problems where a 
stochastic-like resonance can be of interest. One possible 
candidate, that I discuss in this section, is the case of 
Brownian motors. 

The very idea of Brownian motors goes back to Feyn- 
man 



Uj and it can be summarized by saying that we 



want to use the energy of the thermal noise in order to 
make work. This can be done, consistently with the law 
of thermodynamics, if one considers a ratchet potential 
with a periodic time behavior of the temperature, see 
figure (|9|). Going from the top to the bottom of figure 
(|9|), we see that when the temperature is low (cold) a 
particle moving in the ratchet potential is trapped in 
one of the minima. Rising the temperature, the parti- 
cle becomes almost free to move and, because its initial 
condition (minimum) is close the maximum on its right, 
the particle has more probability to go towards the min- 
imum on the right. If the temperature is now decreasing 
the particle has more probability to be trapped in the 
right minimum with respect to be trapped in the mini- 
mum on the left. The overall picture is that of a diffusion 
process with more probability to move on the right than 
on the left. Thus, under the combined action of periodic 
temperature variation and the ratchet potential, we can 
have a non zero drift average velocity of the particle on 
the right direction: we have been able to extract work 
from the temperature field. 

There are many kind of Brownian motors [l2[ , the one 
shown in figure ([9|) is just an example. I will propose 
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FIG. 10: The Brownian motor proposed in section 4. The 
particle moves in a periodic potential and noise is acting on 
it. A perturbing potential — vt) localized in space is 
traveling in the system. 





FIG. 11: The effect of the perturbing potential — vt). 
From the top to the bottom, the particle can flip only in the 
positive x-direction as the potential travels. The flip occurs 
only for a suitable combination of v and the noise amplitude 



here a new kind of motors using the very same idea of 
stochastic resonance to produce a non zero current using 
the noise. Instead of a ratchet potential, let me consider 
a particle which feel a space periodic potential of period 
L. In the limit of overdamped friction, the equation of 
motion of the particle position x is: 



dx 
~dt 



dV 

dx 



V^r]{t) V{x) = Vqcos(2itx/L) (27) 



We now introduce a time dependent forcing F{x — vt) in 
the form of a moving localized perturbation as shown in 
figure ([To]). Actually we can thing of F{^) as the effect 
of a moving potential perturbation, i.e. 



F{x - vt) 



d^{x - vt) 
dx 



(28) 



Without the effect of F, the particle performs a diffusive 
process among the minima of V. The characteristic time 
tl to jump on the right or on the left can be computed 
by using the same method of section 2. The action of the 
localized potential can introduce a new feature which is 
illustrated in figure (pTj) . 

If the speed v of the perturbing potential <l> is properly 
chosen, then the particle has the possibility to jump on 
the right but not on the left when the perturbing poten- 
tial <l> is in a given position. Then the particle remains 
trapped in the new minimum up to the time when the 
potential is now shifted to the next minimum and so on. 
In other words, the particle remains trapped in the mov- 
ing potential $ and it moves "balistically" at the same 
speed V of the perturbing potential. In figure ^T2\i I show 
x{t) with and without $ for a properly chosen value of 
cr. The initial particle position is 300L and the initial 
perturbing potential is localized in x = 0. 

When ^ is close to x = 300, the particle remains 
trapped and it moves with the speed v for a rather long 




500 1000 1500 2000 2500 3000 3500 4000 



FIG. 12: The figure shows x{t) with (straight Une) and with- 
out (dotted line) $(x — vt) for a properly chosen value of a. 
The initial particle position is 300L and the initial perturbing 
potential is localized in x = 0. 



distance. According to our knowledge of stochastic dif- 
ferential equation, at each step the particle is following 
the potential with a probability close to 1 but not exactly 
equal to 1. It follows that the particle is trapped in the 
potential <l> for a distance which is a random variable 
(it is the analogous of the random time tl discussed in 
detail in section 2). 

As for the stochastic resonance, if the noise is to small 
or too large, ^ has no effect on the particle. The same 
is true if the velocity is too small or too large. We can 
check these qualitative statements in the following way. 
We consider an ensemble of N particles and we study the 
average distance {Dx) as a function of noise variance a 
and the non dimensional velocity vtl/L. This is done in 
figures (p!3|) and (fT4j) respectively. 

In both cases, the effect of a resonant like phenomena is 
quite clear. Only for properly chosen values of a and v we 
obtain the maximum of {D^). Some recent experimental 
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FIG. 13: The quantity {Dx) as a function of a 
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non dimensional velocity 

FIG. 14: The quantity {Dx) as a function of v 

results [13'] seem to confirm the numerical analysis done 
in this section. 

5. CONCLUSIONS 



• Stochastic resonance is counter intuitive phe- 
nomenon: it is not trivial that adding noise to a 
system we can enhance the deterministic periodic 
behavior. 

• Stochastic resonance is a robust mechanism ob- 
served in many physical and biological systems. 
The notion of stochastic resonance is now cross 
disciplinary and new applications are found every 
year. 

• We learn a lot in applying stochastic resonance in 
the theory of climatic change. As I mentioned in 
section 2, it is a crucial step to understand that fast 
variables cannot be simply ignored in the study of 
long term climatic change, which is, overall, the 
basic idea introduced by stochastic resonance. 
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In the last 30 years the scientific community did learn 
many different features on complex systems, starting by 
the pioneering works of Lorenz, Ruelle, Mandelbrot and 
others. In many cases, new ideas and tools have been 
introduced in order to "measure" complexity in an ap- 
propriate way. These tools can be used to reveal different 
features of underlying physical or biological mechanisms. 
In some sense, stochastic resonance is also a tool because 
it allows us to understand whether or not non linear ef- 
fects can act in a cooperative way with the complex and 
chaotic behavior of a given system. On the other hand 
stochastic resonance is a mechanism in the full mean- 
ing of the word because it allows to get large effect from 
a small amplitude perturbation. There have been and 
still there are many applications of stochastic resonance 
in problems dealing with the amplification of signal to 
noise ratio, a quite traditional engineering problem. 

In this paper I have reviewed some known and less 
known features of stochastic resonance. Some simple con- 
clusions can be made. 
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